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Aims. We characterize a sample of Gamma-Ray Bursts with low luminosity X-ray afterglows (LLA GRBs), and study their properties. 
Methods. We select a sample consisting of the 12% faintest X-ray afterglows from the total population of long GRBs (lGRBs) with 
known redshift. We study their intrinsic properties (spectral index, decay index, distance, luminosity, isotropic radiated energy and 
peak energy) to assess whether they belong to the same population than the brighter afterglow events. 

Results. We present strong evidences that these events belong to a population of nearby events, different from that of the general 
population of lGRBs. These events are faint during their prompt phase, and include the few possible outliers of the Amati relation. 
Out of 14 GRB-SN associations, 9 are in LLA GRB sample, prompting for caution when using SN templates in observational and 
theoretical models for the general lGRBs population. 
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^ 1. Introduction 

Oh 

i Gamma-ray bursts (GRBs) are the most luminous events in the 
O Universe, with isotropic luminosity between 10 49 - 10 52 erg.s -1 

H ipr~ 


leszaros 


2006| ). GRBs display two comp onents: the prompt 


£2 emission, followed by an afterglow ([Rees & Meszaros||1992 




Meszaros & Rees |1997[|Panaitescu et al. |1998|), b oth observed 


at all wavelengths ( Costa et al. 1997[|Van Paradijs et al. |1997[ 
|Frail et al~||1997| i. In X-rays, the afterglow light curve can be 
described as a steep-flat-steep broken power law (jNousek et al. 
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2006 ). The fir st part (steep decay) has been associated w ith the 
prompt phase ( [Willingale et al. |2007||Zhang et al. |2006| ) while 
the central engine is still active; the rest of the afterglow is due 
to the dynamics of the interaction of the jet with the surrounding 
medium. 

Several studies have been made on GRB samples (e.g. [Me] 
land ri et al. ||2014| ), but in general they do not address specific 
properties. Usually the authors focus on complete samples in or¬ 
der to derive broad properties. 

These properties are then used to define the unknown phys¬ 
ical properties of an archetypal GRB. In this work, we consider 
that the population of long GRBs (hereafter lGRBs) may hide 
various sub-types of GRBs; thus it is important to check for the 
existence of different populations in the sample, and, should this 
happen, how the previous conclusions apply to the whole popu¬ 
lations. This has already been shown with the class of ultra-long 


GRBs (Gendre et al. |2013[ [Zhang et al. |2014[|Boer, Gendre &| 


|Stratta 2014|). 

In the past, several GRBs featuring faint prompt emission 
have been observed: GRB 980 425 ( [Galama et al. |1998l|Kulkarni 


et al. |1998[|Pian et al. |2006[ e.g.), GRB 031203 (e.g.[Malesani 
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-ray bursts: afterglow 


et al. ||2004| [S oderber g et al. |2004[ |Watson et al. |2004), GRB 


060218 (Campana et al. |20 06; Maz zali et al. |200 6 ; Soderberg 


etaI~||2Q06UVii 

et al. ||20nV|Starling et al. ||2011| e.g.). Theoretical work has 


irgili et al. ||2009[ e.g.) and GRB 100316D ([Fan 


been done on the standard model in order to explain these events 


( [Daigne & Mochkovitch|2007]|Barniol-Duran et al. |2015[ e.g.), 

and from their obvious properties (a low luminosity) several au¬ 
thors have pointed out that only very few faint events are de¬ 
tected (e.g. |Imerito et al. |2008| ). All these studies are based on 
the properties of single events (despite the fact that GRB 060218 
and GRB 980425 have very different properties, for instance), 
and so far the global sample of Swift prompt faint events has not 
been studied. Moreover, to our knowledge, despite a different 
approach and properties, no systematic study (and actually no 
individual) has been made by selecting a sample on the basis of 
the faintness of the afterglow. 


This is the purpose of this work, and in the following we 
call GRB displaying dim afterglow (according to our criteria, 
Low Luminosity Afterglow GRBs (LLA GRBs). We explicit the 
criteria to build a consistent sample of LLA GRBs and we derive 
its properties. Then, we use a control sample based on different 
bursts to check whether they form a class different from that of 
normal lGRBs. 


This paper is organized as follows: In section [2] we present 
the LLA GRB sample and we describe how we selected it. In 
section [3] we discus the possible biases and basic properties of 
our sample. In section]?] we discuss our results, before conclud¬ 
ing in section[5] In the following, all errors are quoted at the 90% 
confidence level, and we used a standard flat ACDM model with 
Q m =0.3 and H 0 = 12kms~ l Mpc ~ l . We also used the standard 
notation F oc t~ a v~P. 
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2. Definition of the sample 

We took into account all bursts with a measured redshift ob¬ 
served before 2013, February the 15 th , without consideration of 
the detector triggered by the event and/or observing it. We have 
used the list compiled by GreineQ This leads to a first sample 
of 283 sources which have been observed at X-ray wavelengths, 
including short and long GRBs. As we are interested only in the 
later, we have to exclude sGRBs: to that purpose we used the 
method described in [Siellez et al. [ ( |2014[ ; this method classi¬ 
fies short GRBs all burst with a duration less than 2 seconds in 
the rest frame, with additional criteria on the afterglow) to reject 
them, leaving 254 long bursts in the global sample. 

As the analysis of the bursts that happened prior to 2006 was 
already performed by |Gendre et al. ||2008[ ), we describe here the 
method followed for the Swift bursts only. W e retrieved the XRT 
light curve from the online Swift repository ( [Evans et al. 2009). 

Comparing flux light curves is a complex task, and need a 
careful estimation of the spectral index and the count-to-flux 
conversion factor. The estimation of these two parameters are 
done automatically for the online repository light curve, using 
standard models that may fail for various reasons, or not corre¬ 
spond to our needs (for instance, a spectral index estimated with 
some data taken before the end of the plateau phase, see below). 
We thus cannot use directly the data downloaded from the online 
repository, and needed to estimate independently the spectral in¬ 
dex and the count-to-flux conversion factor. 

For this purpose, we also downloaded the raw data from the 
archives, and applied to them the calibration released in May 
2013 (using the Swift software, distributed as part of the HEA- 
SOFT package, version 6.12). We then extracted a spectrum us¬ 
ing the task xs elect, also part of the HE AS OFT packag^ and 
fit it with a power law model absorbed twice, in the host galaxy 
and in the Milky-Way. The N H value of the Galaxy was set to 
the value given by the Leiden/Argentine/Bonn (LAB) Survey of 
Galactic HI ( |Kalberla et al. |2005| ), while the one for the host was 
let free to vary at the host redshift. Finally we compared this best 
fit model with the one from the automated analysis pipeline: in 
case of inconsistency we recomputed the flux light curve using 
the conversion found from our best fit model. 

Once the flux calibration has been checked and eventually 
corrected, we selected the afterglow part of the light-curve. We 
followed the method of Gendre et al. ( |2008| ), removing from the 
light curves all emission present before the end of the plateau 
phase and all flaring emissions. This net light curve was then 
corrected taking into account the cosmological effects including 
the K-correction. 

We worked on a "flux" light curve, rescaling all light curves 
to a common distance of z = L As stated in IGendre & Boer | 
(2005), this allows a smaller uncertainty on the final light curves. 
One may wonder, now with the precise cosmology parameters 
measured by Planck ( [Hinshaw et al. |[20l3| |Ade et al~||2013| ) 
whether this is really needed; the reason is that the uncertainty is 
introduced by the K-correction (and not by the distance correc¬ 
tion), which is very sensitive to the spectral index: 


Koc(l+zf 


( 1 ) 


As an example, with a redshift of 4 and a precision of 1.0+0.3 
for f3, the uncertainty on K is 5 ± 3, i.e. 60%. Rescaling to z = 1 

1 http ://w w w.mpe. mpg. de/~j cg/grbgen .html 

2 http://www. swift.ac.uk/xrt_curves 

3 http://heasarc.gsfc.nasa.gov/heasoft/ 
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Table 2. Statistical parameters of the cumulative redshift distributions. 


Parameter 

LLA GRBs 

All lGRBs 

mean 

0.55 

2.20 

median 

0.53 

1.98 

standard deviation 

0.38 

1.19 


leads to 2.5±0.7, i.e. an uncertainty of 28% : this method reduces 
the scattering induced by the uncertainties on the measurements, 
allowing for a more precise selection of the sample. 

Being interested on LLA events, we defined two template 
afterglows with a decay index of 1.2 and 1.4 respectively (cor¬ 
responding to the typical values expected with p ~ 2.3 where 
p is the power law index of the accelerated electrons in the 
cases of wind and interstellar media); we set a priori the limit 
at F = 1 x 10“ 13 erg cm -2 s -1 one day after the burst. 

There are two reasons for this choice: first, we are interested 
in the low luminosity part, and thus we chose a flux significantly 
lower than the mean observed flux for the afterglows present in 
Fig. |TJ second, as noted in |Gendre et al. |p008| l, there are bursts 
with a low luminosity that seems to not follow the properties of 
the other groups. These events represent about 10% of the total 
burst population, which turns to a limiting flux of about 10“ 13 
erg cm -2 s -1 one day after the burst. Additionally, the flux cut¬ 
off 10“ 13 erg cm 1 2 3 s -1 at 1 day corresponds roughly to the low¬ 
est afterglow luminosity at one day of the unbiased sample of 
|D’ Avanzo et al. | ( [2012a| ) (there in as seen from Figure 2). 

All bursts with an afterglow light curve entirely below these 
two templates were part of LLA GRBs sample; the others are 
used as a control sample. The result of the selection is displayed 
in Fig. [T] 

The final sample includes 31 events that are listed in Table 

0 representing about 12% of all lGRBs considered here. Table 
displays the GRB name, redshift, galactic and host N#, galac¬ 
tic and host Ay, the afterglow temporal and spectral indexes, the 
isotropic and peak energies, and the T90 duration (the time during 
which 90% of the energy of the prompt is emitted). For those af¬ 
terglows displaying a break after the plateau phase (GRB 060614 
and GRB 120729A), the decay index is indicated pre-break. 


3. Statistical Properties 

3.1. The redshift distribution 


The redshift distributions of the LLA GRB sample is shown in 
Fig. [2] together with the distribution of all lGRBs. 

A simple examination of Fig. [2] shows that LLA GRBs are 
closer than norm al lGRBs whose distribution peaks on averag e 
at z = 2.2 ( e.g. |Jakobsson et al. ||2006[ 


Coward et al. 2013) 


Table [2] displays the main parameters of the two distributions. 
We performed a Kolmogorov-Smirnov test on the two data sets, 
which shows that the probability for the two distributions to be 
based on the same population is 1.1 x 10 -14 , hence rejecting this 
hypothesis. 

We check now whether the difference between the two red¬ 
shift distribution is intrinsic or due to a selection bias. 

Faint events are more difficult to detect than brighter ones. 
Furthermore, the measure of the redshift implies that the op¬ 
tical afterglow is bright enough for spectroscopic observations 
to be performed. As a matter of consequence, LLA GRBs are 
plagued by a detection bias that prevent them to be detected at 
large distance. From their flux, we estimate that the faintest of 
the LLA GRBs present in our sample can be detected up to a 
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Fig. 1 . The light curves of all sources, corrected for distance effects (see text) and rescaled at a common redshift z=l. LLA events are shown with 
blue diamonds, the control sample is shown with red dots. 


distance of z = 1. Because regular lGRBs can be detected up to 
z = 8.2 (Tan vir et al. ||2009] ), we also consider that the sample 
of lGRBs is complete for z < 1, thus removing the detection 
bias. We have recomputed the cumulative redshift distributions 
for this sub-sample (see Fig. [3]). The difference is still large, and 
from a Kolmogorov-Smirnov test the probability that the two 
distributions are drawn from the same population is 9.4 x 10 -4 . 
We thus conclude that the LLA GRB population is different from 
the "classical" lGRBs one. 


3.2. Absorption and Extinction 

For consistency, we first checked that our distribution for the 
Milky Way values of A v ,Goi and N n,Gai (he. the optical extinction 
and X-ray absorption parameters) is consistent with the whole 
sample of normal long GRBs. While the X-ray absorption has 
little effect on our sample since we use the flux in the 2.0-10.0 
keV band, where absorption can be neglected (Morrison et al. 


|| |1983| ), it is well known that the optical extinction can bias a 
distribution (for instance the well known problem of dark bursts, 
e.g. |Jakobsson et al. |2004[ ). 


The optical extinction was calculated using the NASA/IPAC 
extragalactic databas^ 1 2 * 4 ] for the Landolt V band measured by 
[Schlegel et al.~ ( 1998 ) for all bursts but GRB 060904B and GRB 
061110A. These two bursts are seen in projection on the galac¬ 
tic disk where the measures of Schlegel et al. | ( p~998] ) are highly 
variable with the position. For these two events, we rely on the 
most accurate measurements of Sch ady et al. | ( |2012| ) an d|Zafar| 
et al. | ( [2QTT] ) respectively. The results are reported in Fig.|4| One 
can note that the Ay values obtained for LLA GRBs statistically 
are not very larger compare to the one for normal lGRBs (Ay < 
2 for 87% of normal lGRBs, ( [Covino et al. |l2013[ )). In the fol¬ 
lowing, we consider that the gas and dust in the Galaxy have not 
introduced a bias in LLA GRB sample. 

The intrinsic hydrogen column density N h can be linked to 
the host properties ( [Reichart & Price |2002| ), thus we also inves¬ 
tigated on this. The intrinsic host absorptions for the LLA GRBs 


http://ned.ipac.caltech.edu/forms/calculator.html 


are mostly compatible with little or no intrinsic absorption. We 
see that for the sources with a non zero N h,x (Fig. [5]), the absorp¬ 
tion of the host galaxy is on average a factor 10 larger than in the 
Milky Way, as already noted by |Starling et al. | ( |2013| ). At low 
redshift this effect was attributed to the gas in the host galaxy. 


3.3. Decay and spectral index 


The distributions of the temporal decays and spectral index for 
the LLA GRB sample are displ ayed in Fig. [6] We used a ref¬ 
erence sample of bursts listed in |Gendre et al. | ( |2008| and not 
members of the LLA GRB subclass. Note that the decay index 
of GRB 060607A reported in this last article is incorrect and 
is not considered in the comparison. The two samples are very 
similar, as indicated by a K-S test (p = 0.80 and p = 0.08 for 
the decay and the spectral indices respectively). We thus con¬ 
clude that the two samples have similar spectral and temporal 
properties. This can also be seen when considering the closure 
relations (Meszar os et al. |1998l|Sari et al. |1998[|1999a[|Cheva-| 


|lier & Li 2000; Zhang & Meszaros 2008) to investigate the burst 
geometry, the fireball microphysics, its cooling state and the sur¬ 
rounding medium (presented in Fig. [ 7 ]). These are very similar 
to the ones obtain from BeppoSAX, XMM-Newton or Chandra 
( |De Pasquale et al. |2006| |Gendre et al. ||2006| ) for long bursts. 
We note however two peculiar events: 


1. GRB 120729A: The pre-jet break closure relations are re¬ 
jected for this event. We can thus identify the break presents 
at tb = 8.1 ks in the light curve as the jet break, obtaining 
the positions of the specific frequencies and the value of p 
( v m < vxrt < v c , P - 2.8 ±0.2). The pre-jet break decay 
properties are in agreement with this (green point in Fig. [ 7 ]). 

2. GRB 060614: This event would be compatible with another 

jet effect, with p - 2.25 ± 0.05. However, the errors bars are 
large enough to accommodate some non jetted closure rela¬ 

tions. We thus cannot conclude firmly on the jet hypothesis 
for this source based on the closure relations alone. 
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Table 1 . The burst sample and its main characteristics (see text). The spectral and temporal indexes for GRBs before August 2006 are taken from 
|Gendre et al. | ( |2008] ). 


GRB 

z 

Gal Host 

(10 21 cm -2 ) 

Gal 

Ay 

Host 

(mag) 

Afterglow 

Temporal Spectral 

index index 

logT a 

(s) 

F 

^lSO 

(10 52 erg) 

E P ,i 

(keV) 

T 90 

(s) 

Ref. 

GRB 980425 

0.0085 

0.428 


0.071 

1.73 

0.10+0.06 

(0.8) 


(1.3±0.2)xl0 -4 

55+21 

18 

(1), (13) 

GRB 011121 

0.36 

0.951 


0.061 

0.38 

1.3+0.03 

(0.8) 


7.97+2.2 

1060+275 

28 

(1), (14), (15) 

GRB 031203 

0.105 

6.21 


0.117 

0.03 

0.5+0.1 

0.8+0.1 


(8.2±3.5)xl0 -3 

158+51 

40 

(1), (14) 

GRB 050126 

1.29 

0.551 

(0.0) 

0.182 


1 i+ 0 - 6 

05 

0.7+0.7 


[0.4 - 3.5] 

>201 

24.8 

(17) 

GRB 050223 

0.5915 

0.729 

(0.0) 

0.078 

>2 

0.91+0.03 

1.4+0.7 


(8.8±4.4)xl0 -3 

110+55 

22.5 

(2), (18) 

GRB 050525 

0.606 

0.907 

0.38^4 

0.221 

0.36+0.05 

1.4+0.1 

1.1 ±0.4 

3.8 

2.3+0.5 

129+12.9 

8.8 

(5), (19) 

GRB 050801 

1.38 

0.698 

(0.0) 

0.989 

0.3+0.18 

1.25+0.13 

1 84 + 0.56 
1 - oh _0 53 

3.2 

[0.27 - 0.74] 

<145 

19.4 

(5), (17) 

GRB 050826 

0.297 

2.17 


2.398 


1.13+0.04 

1.1 ±0.4 

4.04 

[0.023 - 0.249] 

>37 

35.5 

(17) 

GRB 051006 

1.059 

0.925 

(0.0) 

2.345 


1.69+0.13 

1 c + 0.44 

1 -0 46 

2.77 

[0.9 - 4.3] 

>193 

34.8 

(17) 

GRB 051109B 

0.08 

1.3 

<2 

0.3 


1.1+0.3 

0.7 ±0.4 

3.14 



14.3 


GRB 051117B 

0.481 

0.46 

(0.0) 

0.321 

< 1.4 

1.03+0.5 

(0.8) 


[0.034 - 0.044] 

<136 

9.0 

(11), (17) 

GRB 060218 

0.0331 

1.14 

6+2 

0.437 

0.5+0.3 

1 3 +u 
-0 6 

0.51+0.05 

5.0 

(5.4±0.54)xl0 -3 

4.9+0.49 

-2100 

(1), (20) 

GRB 060505 

0.089 

0.175 

(0.0) 

0.209 

0.63+0.01 

1.91+0.2 

(0.8) 


(3.9±0.9)xl0 -3 

120+12 

~4 

(1), (21) 

GRB 060614 

0.125 

0.313 

0.5+0.4 

0.068 

0.11+0.03 

2.0!^ 

0.8+0.2 

4.64 

0.22+0.09 

55+45 

108.7 

(3), (21) 

GRB 060912A 

0.937 

0.420 

(0.0) 

1.436 

0.5+0.3 

1.01+0.06 

0.6+0.2 

3.3 

[0.80 - 1.42] 

>211 

5.0 

(4), (17) 

GRB 061021 

0.3463 

0.452 

0.6+0.2 

0.185 

<0.10 

0.97+0.05 

1.02+0.06 

3.63 



46.2 

(3) 

GRB 061110A 

0.758 

0.494 

(0.0) 

<0.10 

<0.10 

1.1+0.2 

0.4+0.7 

3.68 

[0.35 - 0.97] 

>145 

40.7 

(3), (17) 

GRB 061210 

0.4095 

0.339 

(0.0) 

0.489 


1.67+0.85 

(0.8) 


[0.10-0.33] 

>105 

85.3 

(17) 

GRB 070419A 

0.97 

0.24 

< 10 

0.081 

<0.8 

0.56+0.0 

(0.8) 


[0.20 - 0.87] 

<69 

115.6 

(5), (17) 

GRB 071112C 

0.823 

0.852 

<5 

0.203 

0 20 +0 05 

1.43+0.05 

0-^4 

3.0 



15 

(4), (17) 

GRB 081007 

0.5295 

0.143 

0 97 +6 - 9 

v y ' - 0.97 

0.196 

0 36 +9 ® 6 

J - 0.04 

1.23+0.05 

n 00+0.88 

u -^_ 0.43 

4.5 

0.18+0.02 

61+15 

10 

(7), (22) 

GRB 090417B 

0.345 

0.14 

22+3 

0.083 

0.8+0.1 

1.44+0.07 

1.3+0.2 

3.54 

[0.17-0.35] 

>70 

>260 

(6), (17) 

GRB 090814A 

0.696 

0.461 

(0.0) 

0.15 

<0.2 

1.0+0.2 

(0.8) 

3.5 

[0.21 - 0.58] 

<114 

80 

(7), (17) 

GRB 100316D 

0.059 

0.82 

(0.0) 

0.088 

2.6 

1.34+0.07 

0.5+0.5 


(6.9±1.7)xl0 -3 

20+10 

292.8 

(8), (23) 

GRB 100418A 

0.6235 

0.584 

(0.0) 

0.623 

0.0 

1.42+0.09 

0.9+0.3 

4.82 

[0.06-0.15] 

<50 

7.0 

(9) (17) 

GRB 101225A 

0.847 

0.928 

(0.0) 

0.311 

0.75 


(0.8) 

4.65 

[0.68 -1.2] 

<98 

1088 

(12), (17) 

GRB 110106B 

0.618 

0.23 

(0.0) 

0.032 


1.35+0.06 

i q 9 + 0.67 
1 - J +_0 32 

4.04 

0.73+0.07 

194+56 

24.8 

(24) 

GRB 120422A 

0.283 

0.372 

(0.0) 

1.241 

0.0 

1.3+0.3 

0.4+0.4 

5.07 

[0.016 - 0.032] 

<72 

5.35 

(10), (25) 

GRB 120714B 

0.3984 

0.187 

(0.0) 

0.077 


1.89+0.02 

(0.8) 


0.08+0.02 

69+43 

159 

(17) 

GRB 120722A 

0.9586 

0.298 

350!™ 

0.555 


1.2+0.4 

1.2+1.2 


[0.51 - 1.22] 

<88 

42.4 

(17) 

GRB 120729A 

0.8 

1.4 

(0.0) 

0.112 


2.8+0.2 

0.8 ±0.3 

3.9 

[0.80 - 2.0 ] 

>160 

71.5 

(17) 


Note: for Ay values: (1) Savaglio et al. (2009); (2) Pellizza et al. (|2006); (3)(Zafar et al. (2011|); (4) Schady et al. 

(2012); (5) 

Kann et al. (2010); (6) Perley et al. (2013); (7) Greiner et al. (2010 

>; (8) Starling et al. (2010); (9) Marshall et al. 

(2011); (10) 

Cano (2013); (11) Michalowski et al. (2012 

); (12) Campana et al. 

(2011); for Ei S0 & E P i values: (13) Pian et al. ( 

1999); (14) 

Ulanov et al. (2005); JT5) Amati et al. (2009 

); (17) ir 

1 this work; XI 

S)[Tabrera et al. (2008); (19) Sakamoto et al. ( 

2011); (20) 

Campana et al. (2006k (21) Amati et al. (2007 ); (22) 

Bissaldi et al. 

(2008); (23) Starling etal. (2011); (24) Bhat (2011); (25) 

E 

[elandri et al. 1 

2012f 




Decay index Spectral index 

Fig. 6. Distribution of the decay index (left) and spectral index (right) for the LLA GRBs sample (blue) compared to the reference sample (red). 


3.4. Prompt phase 

We also investigated the prompt properties of the LLA GRBs. 
For this purpose, as the BAT bandwidth is narrow, we used 
whenever possible the data from Fermi GBM. For events seen 
by Konus-Wind or BeppoSAX, we used previously published 
results. About half of the events have a firm measurement of the 
prompt parameters, the other half presenting upper and lower 
limits. We corrected for the cosmological redshift the values 


of E p , obtaining the intrinsic E Pji values. These values cluster 
broadly within the 40-200 keV range. 

We note however that there is a lack of bright events also 
in the prompt phase. Taking into account the median redshift of 
LLA GRBs, the E pp values, the intrinsic scatter of the Amati 
relation, and the properties of the BAT instrument, one should 
expect to detect bursts up to Ei S0 = 3x 10 53 ergs, at least one order 
of magnitude larger than the brightest measurements listed in 
TableQ] We thus conclude that this effect is an evidence that LLA 
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Fig. 7. X-ray decay index versus spectral index of LLA GRBs. The purple filled circle and square represent GRB120729A and GRB 060614 
respectively. The green dot represent GRB 120729A before the break at t b = 8.1 ks. All closure relations, indicated by the lines, are computed for 
p > 2 in the slow cooling phase. Solid and dash-dotted lines stand for v m < v < v c and v c < v respectively. Blue, red and black lines stand for ISM, 
wind medium, and jet effect respectively. 



Redshift 



Fig. 2. Top: Redshift distribution of LLA GRBs (blue) compared to 
the normal lGRBs (red). Bottom: Cumulative distribution of the same 
samples. 



Fig. 3. Cumulative distribution of the redshift of LLA GRBs (blue) 
compared to all GRBs (red) for redshifts z<l. 


GRBs are intrinsically less energetic, both during the prompt and 
the afterglow phases, compared to normal lGRBs. 


4. Discussion 


4.1. Distance of the sample 


As noted by previous authors (e.g. 

Watson et al. 2004 Guetta 

& Della Valle 2007; Daigne & Moc 

ikovitch 2007; Imerito et al. 

2008), faint GRBs cannot be detected at large distance, and, by 


definition, all LLA GRBs have a faint luminosity afterglow. We 
are thus missing distant LLA GRBs, as one could expect from 
the redshift distribution. 


On the other hand, can we consider that all bursts with a 
faint afterglow are LLA GRBs, and thus that our sample is not 
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A v (mag) 

Fig. 4. Histogram of the optical extinction Ay in the Galaxy for LLA 
GRBs. GRB-SN associations are represented by cyan bars while other 
LLA GRBs are shown with blue bars. 



Fig. 5. The HI column density in the host galaxy, Nhjc for 11 LLA 
GRBs. The blue bars are fitted values while the white ones with red 
arrows are upper limits. 


contaminated by some normal lGRBs? We assume this is not the 
case based on our selection criteria, which allows to discriminate 
regular nearby lGRBs such as GRB 030329 (which is not part of 
our sample, and in fact a normal 1GRB). 

Is this distribution of redshift biased? The volume of the Uni¬ 
verse at low redshift is very small, thus allowing for few events 
to occurs: this could explain the lack of normal lGRBs at redshift 
lower than 0.3. We note however that this argument also apply 
to LLA GRBs, and thus that if the two populations were similar 
in their redshift distribution, we should see the same proportion 
of bursts located between 0 < z < 0.3 and 0.5 < z < 1.0 for LLA 
GRBs and normal lGRBs. As can be seen in Fig. [2] this is not 
the case. We thus conclude that, if our sample is contaminated, 
the proportion of normal lGRBs is not large enough to prevent 
the main properties of this group to be apparent, and that LLA 
GRBs are events closer than normal lGRBs. 


4.2. Geometry and environment of the burst 

Most of the sources can be explained by both a wind envi¬ 
ronment and a constant ISM. As shown below, many of these 
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sources are associated with SNe (see Table |3). Thi s ass ociation 
would point towards a wind environment ( Chevalier et al. 2004). 


However, as shown by Gendre et al. ( 2007| ), the termination 
shock can lie very nearby to the star, and we cannot conclude 
firmly on the surrounding medium. 

One source deserves a more careful study: GRB 120729A. 
This event can be accounted for by the closure relation of a jet. 
There is also a hint of achromaticity, as a break is seen both in 
X-ray and in optical around the same time (Masel lTet al. 12012} 
D’ Avanzo et a l. |2012b| ). This supports the interpretation of a jet 
effect. 

Levinson & Eichler |2005| 


The opening angle is given by 


who extended the work of Sari et a 


radiation efficiency of the prompt phase: 


. |( |1999b| ) to account for the 


0(t b ,E iso ) = 0.161 


( v \ 3/8 

(l+z) ' 


1/8 / % 


-'iso, 52 


1/8 


( 2 ) 


where the standard values for the number density of the medium 
n = 1cm -3 and the radiative efficiency tj 7 = 0.2 are used. We get 
6 = 2.7°. 

From the post jet-break part of the light curve, we derive 
p = 2.8 ± 0.2. This value is not compatible with both the spectral 
and temporal decay indexes (0.74+0.072 and 1.08+0.03 respec¬ 
tively) of the pre-break part of the light curve. Only the spec¬ 
tral index is marginally consistent with this value of p, assuming 
v m < y x < v c and an ISM. The temporal decay is too flat (we 
expect a value of at least 1.5). In order to reconcile all of these 
facts, we need to involve some late time energy injection to flat¬ 
ten the light-curve ( [Hascoet et al. |2014] ). This energy injection 
needs to be present during the pre-break part of the light curve, 
but should stop during the post-break part. We note that the sam¬ 
pling of the X-ray light curve is not good during the jet break 
and allows for some non-simultaneity. 

A similar argument can be drawn for GRB 060614, which 
may also be compatible with a jet, according to the closure rela¬ 
tions. This burst also displays an achromatic break (around 36.6 
ks) in X-ray, optical and UV ( Mang ano et al. 2007| ). However, 
before the jet, this burst features a plateau phase and not a stan¬ 
dard afterglow. This is somewhat unusual for a GRB, and would 
request energy injection fine tuned in order to stop at the moment 
of the jet break. We note in addition that GRB 060614 has been 
proposed to be a short GRB with an extended emission (Zhang 
|et al. ||2007| , making this event clearly odd in our sample. The 
explanation of why the energy injection would stop at the same 
time than the jet break is beyond the scope of this paper. In any 
cases, if we assume the presence of a jet, the corresponding jet 
opening angle is 6.3°. 

Several authors ([Yamazaki et al. [2003} Ramirez-Ruiz et al 


||2005[ |Daigne & Mochkovitch|2007| e.g.) have tried to explain 

the properties of some of these events based on the jet properties 
(aperture angle, viewing angle). From our findings, when we do 
have a measurement of the jet aperture it is similar to the one of 
normal 1GRBs [Ghisellini | ( 201 2[ 6 = 4.7 + 2.3°,), and there is 
no hint of large off-axis viewing. The typical LLA GRB should 
have a jet not very different from the one of normal lGRBs. 


4.3. Microphysics of the fireball 

The spectral and temporal properties, once merged into the clo¬ 
sure relations (see Fig. [ 7 ]), can indicate the position of the cool¬ 
ing frequency and thus give some insight into the parameters of 
the fireball. For LLA GRBs, we have two possibilities: either we 
cannot conclude, or the X-ray band is located below the cooling 
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frequency. In the former case, this is just due to the uncertainties 
of the measurement. In the latter, this is not common: indeed, 
most late GRB afterglows are compatible with the X-ray band 
located above the coo ling frequency (Gend re~et al. ||2006} [De 
|Pasquale et al. |2006 ). We insist here on the fact that this mea¬ 
sure is time dependent, as shown in |Gendre et al. ] (2QQ6), and 
that the comparison need to be done with consistent data. 

In the case of an homogeneous interstellar medium (ISM), 
the formula of the cooling frequency is (|Panaitescu & Kumar 
2000 ): 


v c = 3.7 x 10 14 ^“ 3 1/ V 1 (F+ l)- 2 e“^ 2 2 r; 1 / 2 Hz, 


( 3 ) 


where £53 is the isotropic energy in units of 10 53 ergs, e B ~2 is 
units of 10 -2 , n is the number density of the medium in the units 
of cm ~ 3 , Y is the Compton parameter, e B -2 is the fraction of 
internal energy of magnetic field, Tj is the time expressed in 
days after the burst. 

Instead in the case of a wind medium, the cooling frequency 
reads 


= 3.5 x 10 u eI , 3 2 a; 2 (y +1 r 2 


3/2 t 1/2tt 
e B -2 ^ d ^ Z ’ 


10 ~ A E~l l2 e~^l < 1, 


-4^1/2 -3/2 
"53 t B -2 


may be lower); GRB 031203 may be much softer than measured 
by INTEGRAL/ISGRI as supported by dust echo measured by 
XMM. 

We can indeed see that the outliers are all located in the 
left part of the E p j - Ei so plane. In this part of the diagram, the 
usual gamma-ray instruments are not well suited to measure the 
prompt properties. For instance, the BAT measurements of GRB 
060218 ( [Sakamoto et al. |2006| ) would have make this event more 
similar to GRB 980425, i.e. a clear outliers. It is its simultane¬ 
ous observation by BAT and XRT that makes it compatible with 
the Amati relation. One may thus imagine that this conclusion 
may hold for all outliers. We note however that the prompt phase 
of a GRB usually shows a hard-to-soft spectral variation (e.g. 
Meszaros |2006| ) and that the prompt emission of GRB 060218 


( 4 ) 


, Pji - ^ 1S0 plane changes dramati- 


where A* is the number density in the wind. 

For those bursts where we can conclude on the position of 
the cooling frequency, we can then insert the numbers to have 
an idea of the constraints on the model. We start by assuming 
that the fireball expands in the ISM. The XRT band ranges from 
7.2 x 10 16 Hz to 2.4 x 10 18 Hz respectively. We however assume 
that v c is above 3.7 x 10 18 Hz (i.e. slightly above the X-ray band) 
for simplicity. Equation [3] simplifies to : 


( 5 ) 


when assuming the standard density n = 1cm -3 , the Compton 
parameter Y <tc 1 and considering the observation time of 1 day. 
Taking the lowest £53 measured (in order to have the stringent 
constraint), i.e. 10 -5 (value of £53 for GRB 980425), we obtain: 
€b ~2 > 0 . 1 , which is not really constraining, as typical values of 
e B -2 should be of the order of 1 for lGRBs. 

The situation is similar when assuming a wind medium, for 
which Eq.[4]implies : 


lasted significantly longer than other bursts: it is not sure that a 
measurement consistent with the ones done at the BeppoSAX 
epoch (i.e. time averaged over the complete prompt emission) 
would lead to a similar conclusion. 

On the other hand, GRB 980425, GRB 060505 and 
(marginally) GRB 050826 are not compatible at all with the Am¬ 
ati relation. If we assume that these measurements are correct, 
then the best fit relation in the E nA - Ej 
cally, being far more flatter. 

A flatter Amati relation has been foreseen as early as 2003 
( jYamazaki et al. [20031 ), using GRBs seen off-axis. For complete¬ 
ness, we also note that a similar explanation hold in case of the 
canonball model ( [Dado & Dar[2005| ). Being seen off-axis, these 
events are expected to be less luminous than normal lGRBs even 
during the afterglow (see e.g. D’alessio et al. 2006| . A balance 
need however to be made between this argument and the fact that 
the events are detected: once seen off-axis, the luminosity of the 
burst decrease very fast with the off-axis angle. The simple fact 
that these events are detected means that they are not at a very 
large off-axis angle, but should rather be considered slightly off- 
axis. 

Also, as already indicated, there was no a-priori reason why 
the LLA GRBs have a low Ei S0 (and thus a low prompt lumi¬ 
nosity). If it is known that the prompt and afterglow luminosi- 
ties are linked together (|De Pasquale et al. |2006} |Gehrels et al. 


( 6 ) 


when assuming a standard density and a Compton parameter, 
and considering the observation time. Using the same method, 
but this time using the largest value of £53 (again in order to have 


2008J |Evans et al. ||20Q9| ), this result holds for normal lGRBs 

and should not considered as obvious for LLA GRBs. We can 
confirm it also holds for this class of events. 

4 . 5 . LLA GRBs and supernovae 

Nine LLA GRBs are firmly associated to SNe by spectral and 


the stringent constraint), we obtain e B -2 > 0.0012. Again, the 

gether with the other solid associations 

0 

magnetic energy of the fireball seems not to explain the unusual 

In our sample, GRB 060505 ( Haisl 

lip et al. 2006) and GRB 

position of the cooling frequency. 

060614 (gynbo et al. 2006; Della Va 

lie 2007; Gal-Yam et al. 


mon position of the cooling frequency for LLA GRBs compared 
to normal lGRBs is due to the small energy of the fireball rather 
than the magnetic energy of the fireball. 

4 . 4 . Prompt properties of LLA GRBs 

In Fig. [ 8 ] we clearly see that all the outliers of the Amati rela¬ 
tion belong to the LLA GRB sample. Several explanations have 
been proposed to explain these events (see |Amati ||2006| |Am 


distance, any SNe would have been detected, and thus these non 
associations are significant. This may strengthen the conclusion 
that these two events are two short bursts, or at least two events 
not related to a normal colapsar. From our criteria, these events 
are not short, and thus clearly belong to the LLA GRB subclass. 


|ati et al. |2007[ and reference therein for details): GRB 060505 
may be a short GRB (as its location in the E p j - Ei S0 plane may 
suggest); the E p j value of GRB 061021 refers to the first hard 
pulse, while a soft tail is present in this burst (so the true E p j 


5 While we consider in the following only positive associations, we 
note that two other sources (GRB 070419A, GRB 100418A) might be 
associated to SNe. GRB 070419A displa ys a faint bump in its light 
curve similar to the one of GRB 980425 ( [Hill et aL"]|2007[ ); A bright 
host galaxy may prevent to observe the signature of a faint SN associ¬ 
ated with GRB 100418A, which can be fainter than r magnitude -17.2, 
comparable to the magnitude of the most faintest Ic SNe (INiino et al. | 

m2) . 
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E iso (erg) E iso (erg) 

Fig. 8. Location in the E A ;-E iso plane of LLA GRBs sample. Left, comparison of both short and long GRBs. Right, compared to short events, with 
firm measurements and non-compatible lower limits on E Pti only. 


We then can conclude that at least some of these events can be 
explained by a different kind of progenitor compared to normal 
lGRBs. 

We also note that a large fraction of the GRB-SNe asso¬ 
ciations (64%) belong to the LLA GRBs sample. The posi¬ 
tive associations include several well known events, such as 
GRB 980425/SN1998bw, GRB 031203/SN20031w, and GRB 
060218/SN2006aj. We stress that this can lead to some prob¬ 
lems, as GRB 980425/SN1998bw is commonly used as template 
for light curves and spectra when looking for a SNe within the 
dataset of a given GRB. As discussed above LLA GRBs progen¬ 
itors may differ from normal 1GRB ones, which also applies to 
the associated supernovae. It is thus more accurate and safe to 
use as template GRB 030329/SN2003dh for normal lGRBs. 


5. Conclusions 

As already noted all papers published so far are based on the ob¬ 
served properties of singular faint events, or in some rare cases, 
two GRBs. GRB 980425-like bursts are faint in their prompt 
phase and cannot be detected at large distance; in addition they 
are not compatible with the Amati relation linking E iso to E p . 
However, this does not extend to all burst with prompt low lu¬ 
minosity that are outliers of the Amati relation. As an example, 
GRB 060218 does follow this relation, despite being faint in its 
prompt phase. In order to have the global answer, one needs to 
perform the study on a statistically significant sample of events, 
and this has not been done so far to our knowledge. 

Here we have based our study and selected the sample on 
the afterglow properties, after the plateau phase, alone. The fact 
that this made incidentally entering some GRBs with a low- 
luminosity prompt was therefore not pre-supposed. This ap¬ 
proach allows us to build a large sample, allowing for statistical 
studies and deriving global conclusions. 

We find that LLA GRBs are on average closer than normal 
lGRBs. Both Ay and N# of LLAs are similar to those of nor¬ 
mal lGRBs. Most LLA GRBs are consistent with the closure 
relations expected by the fireball model. The few outliers can be 
accounted for by an early jet break. We show evidences that the 
events in our LLA sample are also intrinsically fainter during 
their prompt phase, reinforcing the evidence for a different pop¬ 
ulation. Actually, some events do not follow at all the E p j - Ei S0 
relation. 
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We have shown that in order to explain all of these properties, 
we can involve either the geometry of the bursts or a different 
kind of progenitors. In the former hypothesis, the bursts would 
be seen slightly off-axis in order to explain the low energy bud¬ 
get observed in these events. In the latter, one could imagine that 
the progenitor provides less mater for the accretion, thus dimin¬ 
ishing the energetic budget at start. 

The LLA GRBs sample includes a significant fraction of 
the supernovae associated with GRBs (64%), including the best 
known associations. This means that the conclusion drawn on 
the general GRB-SN association is based on a sub-sample of 
the low-luminosity population (the LLA GRBs) that might not 
be representative. We stress the need to confirm this point 
and the previous work on GRB-SN associations using different 
spectral and light curve templates, for instance those of GRB 
030329/SN2003dh. 
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Table 3. The table of GRB-SN associations. The L iso values are calculated from the E iso values which are given in the Table]!]. 


GRB 

name 

redshift 

SN 

identification 

SN 

name 

SN 

type 

Host 

type 

LLA 

GRB 

Li so of GRB 
(10 49 erg s 1 ) 

Ref. 
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0.0085 

spectral 
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BL-lc 
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(SbcD) 
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(1) 
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0.36 
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SN2001ke 

Iln 

N/A 

yes 
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(2) 

GRB 021211 

1.01 

spectral 

SN20021t 

~Ic 

N/A 

no 

969 

(3), (4), (5) 

GRB 030329 

0.168 

spectral 

SN2003dh 

BL-lc 

N/A 

no 

75 

(6), (7), (8) 

GRB 031203 

0.105 
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SN20031w 

BL-lc 

Irr 

Wolf-Rayet 

yes 

0.56 

(9), (10) 

GRB 050525 
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spectral 

SN2005nc 
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417 

(11), (12) 
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0.0331 

spectral 

SN2006aj 

BL-Ib/c 
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bump 
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Ic 
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30 

(16), (17) 

GRB 091127 

0.49 

bump 

SN2009nz 

BL-lc 

N/A 

no 
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(18), (19) 
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spectral 

SN2010bh 

BL-lc 

Spiral blue 
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(20), (21) 

GRB 101219B 

0.55 

spectral 

SN2010ma 

Ic 

N/A 

no 

29 

(22), (23) 
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0.283 

spectral 

SN2012bz 

Ib/c 

N/A 

yes 

0.44 

(24), (25) 

GRB 120714B 

0.3984 

spectral 

SN2012eb 

I 

no identification 

yes 

0.7 

(26), (27) 

GRB 130215 

0.597 

spectral 

SN2013ez 

Ic 

N/A 

no 

75 

(28), (29) 


Note: for GRB-SN associations: (1) Galama et al. (j 1998); (2)|Bloom et al. (2002); (3) Crew et al. ( 

2002); (4) Della Valle et al. 
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